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Abstract 


Marengo-Mogollón, H., Aldama-Rodríguez, A., $: Romero- 
Castro, I. Composite Roughness in Hydraulic Models. Water 
Technology and Sciences (in Spanish), 7(6), 159-178. 


An experimental survey study was carried out on hydraulic 
models in tunnels working as full pipe on composite 
roughness; the experimental analysis was made with 
different friction coefficients (that belong to different 
materials) and compared with values calculated from a 
proposed theoretical formulation. The model was made for 
the diversion work of El Cajón hydroelectric project, these 
criteria also was applied to the Grijalva Tunnels under 
operation in Mexico and they are also being used for La 
Yesca hydroelectric project presently under construction. 


Keywords: Acrylic-plastic-carpet, coefficient of resistance, 
composite roughness, hydraulic models. 


Introduction 


With an analysis by overtopping of the diversion 
works of Aguamilpa Dam, in Mexico, on January 
1992, Marengo-Mogollón (2006) concludes that: 
“the overtopping event would have been avoided if 
a hydraulic concrete lining would have been built at 
the floor, and shotcrete at the walls and vault of the 
tunnels with vault type (16 x 16 m) section (composite 
roughness concept) even though the peak inflow rate 
exceeded by 50% the original design value”. This 
paper shows the analysis made in a hydraulic 
model with composite roughness that simulates 
permanent flow in tunnels working as full pipe. 
While the flow in prototype is mainly governed 
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Marengo-Mogollón, H., Aldama-Rodríguez, A., € Romero-Castro, 
LI. (noviembre-diciembre, 2016). Rugosidad compuesta en modelos 
hidráulicos. Tecnología y Ciencias del Agua, 7(6), 159-178. 


En un aparato experimental se llevó a cabo un modelo en túneles 
hidráulicos trabajando a tubo lleno con rugosidad compuesta; 
el análisis experimental fue con distintos coeficientes de fricción 
(que pertenecen a distintos materiales). Se comparan con valores 
calculados a partir de formulaciones teóricas. El modelo se hizo para 
la obra de desvío del proyecto hidroeléctrico El Cajón, y se aplicó 
también para el diseño de los túneles de emergencia del Grijalva, 
ambos proyectos en operación en México. 


Palabras clave: acrílico-plástico, coeficiente resistencia, rugosidad 
compuesta, modelos hidráulicos. 
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by gravity, in the model it is also influenced by 
viscosity nevertheless, working within a limited 
range of variation of the Reynolds number, and 
adopting Froude's similitude, has made pos- 
sible to ignore such influence in all the analyzed 
models. The main objective of this paper is to 
validate the theoretical analysis made by Elfman 
(1993), Marengo-Mogollón (2005) and to prove 
the best criteria within five formulas that permit 
the estimation of the roughness coefficient. This 
paper is organized as follows: first, a brief de- 
scription of the experimental apparatus is made, 
then is showed a brief hydraulics review of flow 
resistance equations, the hydraulic theoretical 
development in order to evaluate the composite 
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roughness is presented, and the paper concludes 
by comparing the hydraulic experimental 
results of the Colebrook equation with the 13 
criteria in the transition zone. 


Experimental Apparatus 


The experimental apparatus is shown in figure 1. 
The test section is 0,133 x 0,133 m, and has vari- 


Figure 1. Experimental apparatus. 


able slope with a 9m length. The inlet geometries 
tested in the model and built in prototype, are 
rounded. The testes were made with each mate- 
rial acrylic, sandpaper, plastic and carpet (figure 
2), in order to know the main hydraulic proper- 
ties in each one of them. It was then tested a tun- 
nel of compound roughness that was obtained 
when it was used the acrylic in the bottom, sand- 
paper, plastic and carpet in the walls and vault 
respectively in each one of them. In the analysis, 
it was considered that Froude's similitude law is 
applicable, considering that, the friction factor is 
independent of Reynolds' number. 


Resistance and Roughness Coefficients 


Leopardi (2004) stated that the tractive force t 
produced by the current flow between two sec- 


tions is proportional to energy gradient, l = = 
and is expressed as: 


T=yRI (1) 
At turbulent flow state, t depends on density 


p (p =7/8) mean velocity V, hydraulic radius R 
and the roughness K, that is: 


Acrylic 


«ue 
a 
Ml 
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Figure 2. Material used. 
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1=f(p,V,R,K) (2) 


Which, using Buckingham's theorem, can be 
expresses as: 


A 72 
=*pV 3 
ge (3) 


Substituting (3) in (1) makes therefore pos- 
sible to express the energy gradient as: 


E NE (4) 

gR 

The head loss between two sections can be 

calculated by solving numerically the momen- 

tum equation (dh /ds) + A = 0, where H is the 

trinomial of Bernoulli. Therefore, between two 
sections the following relation holds: 


2 2 
2 +h +0 L=2+1,+0,2+A (5) 
8 8 


where the index 1 and 2 specify the considered 
section. 

In order to determine the hydraulic gradient 
along the tunnels, piezometers were installed 
in the sections 6D and 28D of the models 
(L, ,= 3 127 m). According to the entrance re- 
cords upstream from section, 6D there are strong 
effects of contraction and in all the experiments, 
downstream from section 28D there is a clear ef- 
fect of air entrance; from this section, most of the 
hypotheses made for the hydraulic functioning 
as full pipe are no longer valid. 

From equation (3): 


81 _8gRI 


a 6 
py? av? e) 


Using the above relations makes it actually 
to calculate 1, 1, A, for different values of dis- 
charge in each model. 

Also is calculated the absolute roughness 
with Colebrook-White criteria; for R, > 25,000 
the Colebrook-White (Yen, 2002) relation in the 
transition zone is often used, 


K K 


(el + 3 


K,R —4R 


-K;, log (7) 


¿La 
JA 


For full circular pipe Colebrook (1939), 


K,=2.00, K,= 14.83, K,= 2.52, where R,= A 


v 
(V = cross sectional average velocity, and 
R = hydraulic radius). In a vault type section 
(figure 3). 


2 
APA 0, 
P.— b(4+x) 
Anda O 
-  Zy (e) 


Then (Eq. (7)) for vault type sections stays: 


1 2K 2.52 
EAS -2 log Ed (8) 
NN 14.83b 2Vb NN 
V 


Eq. (8) will be used in the analysis. 


Hydraulic Models 


Models with One Material 


Table 1 summarizes the hydraulic parameters 
calculated for acrylic tunnel model (slopes A1- 
0.0007, A2-0.001, A3-0.004), for each one, itis re- 
ported the slope, the flow rate Q (that was tested 
in the range 0.01 < Q < 0.025 m?/s), the mean 
velocity V = Q/A, the Reynolds number Re in 
the range 8 x 10*< R <2.4 x 10” (with kinematic 
viscosity of water v = 0.000001 m?/s), then is 
calculated for each case 1, 1, A, and the absolute 
roughness K, for each criteria (K, , Nikuardse, 
K,,, Halana, K, .,, Churchill and K, Swamee, see 
table A.1), that is used like bottom roughness in 
the analysis. 

Table 2 shows the hydraulic parameters for 
sand paper (slopes S1 - 0.001, S2 - 0.004, 
S3 - 0.008), table 3 for plastic (slopes P1 - 0.001, 
P2 - 0.004, P3 - 0.008) and table 4 for carpet 
(slopes C1 - 0.001, C2 - 0.004, C3 - 0.008). 

For each model, they are reported the same 
parameters of table land the absolute roughness 
measured with Colebrook-White formula; E 
means absolute roughness measured in each 
material (K 
carpet). 


mC,S sandpaper, 0 plastic, E 
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Wall 


2 
A=2b? ¿Lam =0.01579 m? 
b P=4b+nb=b(4+T)=0.4749 m 


P, =2b=0.133 m 


P.,=2b+mb=b(2+m)=0.34192 m 


Pp b PD 25m 
2b(4+1) 2 
Bottom 
Figure 3. Geometry properties of vault type section. 
Table 1. Model tests. Acrylic. 
a td o es $ E dea e A en de fa 

A11 0.0192 1.216 162 077 0.0099 0.3278 0.01740 | 0.000080 | 0.000027 | 0.000063 | 0.000022 

A12 0.0214 1.359 180 648 0.0126 0.4183 0.01787 | 0.000090 | 0.000044 | 0.000075 | 0.000037 

A1 0.0007 A13 0.0228 1.444 192 467 0.0143 0.4741 0.01785 | 0.000089 | 0.000046 | 0.000075 | 0.000040 
Al4 0.0245 uba 206 817 0.0157 0.5208 0.01698 | 0.000072 | 0.000030 | 0.000058 | 0.000025 

A15 0.0258 1.634 217 791 0.0181 0.6007 0.01766 | 0.000085 | 0.000047 | 0.000072 | 0.000041 

A21 0.0196 1.241 165 454 0.0104 0.3459 0.01762 | 0.000084 | 0.000034 | 0.000068 | 0.000028 

A22 0.0214 1.355 180 648 0.0119 0.3972 0.01697 | 0.000072 | 0.000024 | 0.000057 | 0.000019 

A2 0.001 A23 0.0229 1.450 193 311 0.0140 0.4650 0.01735 | 0.000079 | 0.000035 | 0.000065 | 0.000029 
A24 0.0245 1.552 206 817 0.0153 0.5103 0.01664 | 0.000065 | 0.000024 | 0.000052 | 0.000019 

A25 0.0257 1.628 216 947 0.0177 0.5887 0.01744 | 0.000081 | 0.000042 | 0.000068 | 0.000036 

A31 0.0228 1.444 192 467 0.0139 0.4632 0.01744 | 0.000081 | 0.000037 | 0.000066 | 0.000031 

A32 0.0242 1.533 204 285 0.0160 0.5325 0.01779 | 0.000088 | 0.000047 | 0.000075 | 0.000041 

dd A33 0.0261 1,653 220 324 0.0182 0.6064 0.01742 | 0.000080 | 0.000042 | 0.000068 | 0.000036 
A34 0.0275 1.741 232 142 0.0201 0.6667 0.01725 | 0.000077 | 0.000041 | 0.000065 | 0.000035 

0.000080 | 0.000037 | 0.000066 | 0.000031 
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Models with composite roughness 


Table 5 shows the hydraulic properties of 
acrylic-sand paper measured in two roughness 
model (slopes AS1 - 0.001, AS2 - 0.004, AS - 
0.008), table 6, acrylic-plastic (slopes AP1 - 0.001, 
AP2 - 0.004, AP3 - 0.008), and table 7, acrylic 
carpet slopes (AC1 - 0.001, AC2 - 0.004, AC3 - 
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0.008). For each model, it's calculated like before 
L í, A and the measured absolute roughness 
Kc; calculated with Colebrook-White criteria 
cas 
for acrylic plastic and K,,. ,. 
general, this calculus is taken like the “true 


is measured for acrylic-sandpaper, K,,. ¿p 


acrylic-carpet). In 


value” in the transition zone and itis used for 
comparison in all over analysis. 
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Table 2. Model tests. Sandpaper. 


d a O a se z (egjmr) a de 
S1 0.01700 1.077 143 183 0.0629 2.09213 0.03295 0.00082 
S12 0.01900 1.203 160 028 0.0683 2.27077 0.03204 0.00075 
E S13 0.02000 1.267 168 450 0.0743 2.47121 0.03379 0.00090 
S14 0.02100 1.330 176 873 0.0793 2.63773 0.03387 0.00091 
S15 0.02300 1.457 193 718 0.0848 2.82051 0.03377 0.00091 
S21 0.01600 1.013 134 760 0.0579 1.92478 0.03460 0.00096 
522 0.01800 1.140 151 605 0.0625 2.07912 0.03417 0.00093 
S2 523 0.01900 1.203 160 028 0.0678 2.25560 0.03467 0.00098 
0,004 S24 0.02100 1.330 176 873 0.0734 2.43989 0.03440 0.00096 
825 0.02200 1.393 185 295 0.0793 2.63684 0.03452 0.00097 
526 0.02400 1.520 202 140 0.0854 2.84097 0.03428 0.00095 
S31 0.01800 1.140 151 605 0.0622 2.06781 0.03656 0.00116 
S32 0.02000 1.267 168 450 0.0682 2.26757 0.03455 0.00097 
po S33 0.02100 1.330 176 873 0.0740 2.46154 0.03622 0.00113 
S34 0.02300 1.457 193 718 0.0794 2.64104 0.03647 0.00116 
S35 0.02400 1.520 202 140 0.0858 2.85435 0.03657 0.00117 
Table 3. Model tests. Plastic. 
, E celo cu de a (gl) ñ dp 
P11 0.01500 0.950 126 338 0.01342 0.44636 0.038822 0.00138 S 
P12 0.01700 1.077 143 183 0.01719 0.57149 0.038698 0.00137 2 
Re " P13 0.01900 1.203 160 028 0.02136 0.71020 0.038499 0.00136 E 
P14 0.02000 1.267 168 450 0.02376 0.79011 0.038655 0.00137 2 
P15 0.02100 1.330 176 873 0.02585 0.85946 0.038139 0.00132 E 
p21 0.01300 0.823 109 493 0.01026 0.34103 0.039490 0.00145 E 
P22 0.01500 0.950 126 338 0.01384 0.46014 0.040021 0.00151 E 
P2 P23 0.01700 1.077 143 183 0.01765 0.58679 0.039734 0.00149 E 
an P24 0.01800 1.140 151 605 0.01946 0.64709 0.039084 0.00142 3 
P25 0.01900 1.203 160 028 0.02214 0.73605 0.039900 0.00151 pe 
P26 0.02100 1.330 176 873 0.02595 0.86269 0.038282 0.00134 E 
P27 0.02200 1.393 185 295 0.02885 0.95918 0.038782 0.00139 > 
P31 0.01800 1.140 151 605 0.01906 0.63388 0.038286 0.00133 a 
P32 0.01900 1.203 160 028 0.02197 0.73037 0.039592 0.00147 = 
od P33 0.01900 1.203 160 028 0.02174 0.72283 0.039184 0.00143 E 
P34 0.02000 1.267 168 450 0.02401 0.79821 0.039051 0.00142 E 
P35 0.02200 1.393 185 295 0.03063 1.01834 0.041174 0.00166 S 
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Table 4. Model tests. Carpet. 


úl ed crío colo dl E (egin) ll a 
cía 0.00970 0.614 81 698 0.013243 0.44032 0.091582 0.02191 
Cc12 0.01080 0.684 90 963 0.017006 0.56546 0.094871 0.02344 
cl C13 0.01190 0.754 100 228 0.020407 0.67854 0.093769 0.02294 
Di C14 0.01260 0.798 106 124 0.023445 0.77955 0.096091 0.02402 
C15 0.01350 0.855 113 704 0.026529 0.88208 0.094715 0.02339 
C16 0.01430 0.906 120 442 0.030292 1.00721 0.096389 0.02417 
C21 0.00900 0.570 75 803 0.012071 0.40137 0.096970 0.02440 
C22 0.01000 0.633 84 225 0.015517 0.51595 0.100969 0.02629 
C2 C23 0.01100 0.697 92 648 0.018555 0.61697 0.099783 0.02574 
0.004 C24 0.01200 0.760 101 070 0.021956 0.73004 0.099212 0.02548 
C25 0.01300 0.823 109 493 0.025856 0.85970 0.099550 0.02565 
C26 0.01400 0.887 117 915 0.028395 0.94413 0.094266 0.02318 
(27 0.01500 0.950 126 338 0.031977 1.06324 0.092476 0.02236 
C31 0.01000 0.633 84 225 0.015300 0.50874 0.099557 0.02563 
C32 0.01150 0.728 96 859 0.019064 0.63388 0.093797 0.02295 
C33 0.01200 0.760 101 070 0.021785 0.72434 0.098437 0.02512 
pe C34 0.01350 0.855 113 704 0.025503 0.84797 0.091052 0.02170 
C35 0.01400 0.887 117 915 0.028586 0.95049 0.094901 0.02348 
C36 0.01450 0.918 122 127 0.031624 1.05150 0.097871 0.02486 
(37 0.01500 0.950 126 338 0.035070 1.16609 0.101421 0.02654 
Table 5. Model tests. Acrylic-Sandpaper. 
E 
eN 
a : Code 00 e de a NoE sy e : 
z AS11 0.01700 1.077 143 183 0.01206 0.40113 0.02716 0.000404 
a 
Y AS12 0.01900 1.203 160 028 0.01515 0.50365 0.02730 0.000418 
E AS1 0.001 AS13 0.02000 1.267 168 450 0.01682 0.55943 0.02737 0.000424 
E AS14 0.02200 1.393 185 295 0.02023 0.67251 0.02719 0.000417 
E AS15 0.02300 1.457 193 718 0.02345 0.77955 0.02884 0.000523 
E AS21 0.01800 1.140 151 605 0.01466 0.48731 0.02943 0.000553 
E AS22 0.02000 1.267 168 450 0.01783 0.59284 0.02900 0.000528 
E AS2 0.004 AS23 0.02100 1.330 176 873 0.01937 0.64411 0.02858 0.000502 
= AS24 0.02200 1.393 185 295 0.02214 0.73607 0.02976 0.000584 
E AS25 0.02400 1.520 202 140 0.02468 0.82050 0.02788 0.000462 
g AS26 0.01800 1.140 151 605 0.01666 0.55397 0.03346 0.000866 
sl AS31 0.02000 1.267 168 450 0.01888 0.62785 0.03072 0.000650 
E qn AS32 0.02100 1.330 176 873 0.02178 0.72434 0.03214 0.000763 
3 AS33 0.02300 1.457 193 718 0.02682 0.89169 0.03299 0.000836 
g AS34 0.02400 1.520 202 140 0.02863 0.95200 0.03234 0.000784 
E 
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Table 6. Model tests. Acrylic-Plastic. 


d pad 2 a eS e aa dl E 
AP11 0.01600 1.013 134 760 0.01299 0.43186 0.033013 | 0.000822 
AP12 0.01700 1.077 143 183 0.01565 0.52023 0.035227 | 0.001023 
AP20.001 |  AP13 0.01900 1.203 160 028 0.01891 0.62879 0.034086 | 0.000923 
AP14 0.02000 1.267 168 450 0.02177 0.72377 0.035409 | 0.001047 
AP15 0.02100 1.330 176 873 0.02431 0.80820 0.035864 | 0.001092 
AP21 0.01700 1.077 143 183 0.01579 0.52500 0.035550 | 0.001053 
AP22 0.01800 1.140 151 605 0.01837 0.61094 0.036900 | 0.001188 
AP3 0.004 |  AP23 0.02000 1.267 168 450 0.02109 0.70140 0.034315 | 0.000946 
AP24 0.02100 1.330 176 873 0.02345 0.77980 0.034604 | 0.000974 
AP25 0.02300 1.457 193 718 0.03121 1.03761 0.038385 | 0.001350 
AP26 0.01800 1.140 151 605 0.01743 0.57964 0.035010 | 0.001005 
AP31 0.02000 1.267 168 450 0.0211 0.70176 0.034333 | 0.000948 
AP40.008 |  AP32 0.02100 1.330 176 873 0.02419 0.80429 0.035690 | 0.001075 
AP33 0.02200 1.393 185 295 0.02673 0.88872 0.035933 | 0.001100 
AP34 0.02400 1.520 202 140 0.03174 1.05547 0.035859 | 0.001096 
Table 7. Model tests. Acrylic-Carpet 
x Eres cd cs $3 z (g/m di dEl 
ACI1 0.0115 0.728 96 859 0.01415 0.47048 0.06962 0.00623 > 
ACI2 0.0124 0.785 104 439 0.01633 0.54285 0.06909 0.00613 a 
AC 0.001 ACI3 0.0136 0.861 114 546 0.02018 0.67100 0.07099 0.00651 El 
ACI4 0.0145 0.918 122 127 0.02313 0.76900 0.07158 0.00663 S 
ACI5 0.0154 0.975 129 707 0.02603 0.86549 0.07142 0.00660 2 
AC?21 0.0163 1.032 137 287 0.02893 0.96198 0.07086 0.00649 E 
AC22 0.0110 0.697 92 648 0.01293 0.43001 0.06955 0.00621 3 
AC23 0.0120 0.760 101 070 0.01602 0.53254 0.07237 0.00678 E 
cea AC24 0.0130 0.823 109 493 0.01937 0.6441 0.07458 0.00724 E 
AC25 0.0140 0.887 117 915 0.02146 0.71346 0.07123 0.00656 S 
AC26 0.0150 0.950 126 338 0.02450 0.81447 0.07084 0.00648 E 
AC31 0.0130 0.823 109 493 0.01870 0.62181 0.07200 0.00671 E 
AC32 0.0140 0.887 117 915 0.02178 0.72434 0.07232 0.00678 E 
AC33 0.0150 0.950 126 338 0.02519 0.83741 0.07283 0.00689 El 
AC3 0.008 3 
AC34 0.0160 1.013 134 760 0.02859 0.95049 0.07266 0.00685 S 
AC35 0.0170 1.077 143 183 0.03181 1.05754 0.07161 0.00664 E 
AC36 0.0180 1.140 151 605 0.03471 1.15403 0.06970 0.00627 el 
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Composite Roughness Theoretical 
Development 


In an analysis of Head Loses in Tunnels with 
two roughness (Elfman, 1993) used Niku- 
radse Equations; in this article, are used other 
roughness criteria like Haland, Churchill, and 
Swamee; also are considered the (16) Eqs. Cri- 
teria proposed by Yen (2002) for the analysis of 
acrylic-plastic and acrylic-carpet (the Reynolds 
number are very high and turbulence is fully 
developed). 

A tunnel that works with composite rough- 
ness is shown in figure 4; has a length L with a 
cross section that is divided into the correspond- 
ing areas A, and A, delimited by the perimeters 
P, and P 
mid part of tunnel “c” and the contours of equal 


.,, has the point of maximum velocity at 
velocity intersect themselves at right angles. 

It is considered that there are two types of 
roughness K, (absolute roughness size grain 
material) at the bottom and K, in the walls and 
vault. Total shear force F acting along the tunnel 
surfaces is equal to the sum of shear force at the 
bottom (F, = LP,t,) and shear force at walls and 
vault (F,, = LP 1, ); like F=F, + F,: 


> LP1=LP,1, +LP,r,, (9) 
Considering 


1=pV?%1/8; 1, =pV?a,/8; 1, =pV?h,/8 (10) 


wo w 


PV?A=P,V¿h, +P,V?A (11) 


www 


Maximum velocity V 


max 


(point “c”, figure4) 
is equal to maximum velocity V,., of the tunnel 
bottom and the maximum velocity at walls and 
vault V 


max w" 


vd. (12) 


max maxb max10 


From Nikuradse (1933) analysis: 


max 


1/2 
V -v+37s[ (13) 
p 


For the experimental analysis, form Eq. (12) 
and (13) in Eq. (11): 
P P, 


1/2 2 an a 7+ 1/2 = 7 (14) 
(+1,326) (4, +1,326)  (2,17+1,326) 


Figure 4. Cross section of the tunnel and contours of equal velocity. 
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Head losses at the bottom are equal to those 
at the walls and vault: 


2 Di 

Ah, ss Ah, = L,V, A, a En, Vaz Nay (15) 

22R, 2eR, 
. Vide an Vi, ACE pe PV oh, (16) 

y A, As : A, A, 

P, E, 

Like: V, = Vias and V, = mao 
(1+1.3262,”) (1+1.3261.,”) 
=> E, = E (17) 


A,(%;+1,326)" A,(»;”+1,326) 

There are two scenarios of analysis: 

Scenario 1. If the total roughness factor 2 is 
experimentally known, and also the material of 
the tunnel bottom K,, form Eq. (17): 


In Eq. (11): 
AAA > de (19) 
(1 +1,326)  (a,'?+1,326) b 


Since A=A, +A: 


ad ] E, (20) 


A(n"”+1,326) A,(A7”+1326) 


Tables 8, 9, 10 and 11 summarizes the calcu- 
lus of the variables showed in this first scenario 
As, A, My My K.,;; K,,, means the absolute rough- 
ness in walls and vault (í, is N for Nikuradse, H 
for Haland, Ch for Churchill and S for Swamee). 

Like first comparison, K,,, calculated with 
Nikuradse, Haland, Churchill and Swamee 
criteria is compared with measured sandpaper 
model (figure 5). In general the data fits very 


10,17 


well the behavior of measured values (table 3), 


At F and calculated values (tables 8, 9, 10, and 11). 


A A A 
A,(n;?+1,326)' (2,2 +1,326)' pre 


Table 8. Calculus of Acrylic-sandpaper (Nikuradse criteria), Scenario l. 


> ar CA s A 8 da dl le das Pe 
AS11 0.017 0.027162 1.077 143183 0.00302 0.0190 0.03052 0.000762 
AS12 0.019 0.027302 1.203 160028 0.00301 0.0191 0.03071 0.000778 
0.001 AS13 0.020 0.027369 1.267 168450 0.00300 0.0191 0.03081 0.000787 


AS14 0.022 0.027191 1.393 185295 0.00302 0.0190 0.03056 | 0.000765 


AS15 0.023 0.028838 1.457 193718 0.00288 0.0193 0.03283 | 0.000976 


AS21 0.018  |0.0294330 1.140 151605 0.00283 
AS22 0.020 [| 0.0290041 1.267 168450 0.00286 


0.0193 0.03366 | 0.001059 


0.0193 0.03306 | 0.000998 


0.004 AS23 0.021 0.028582 1.330 176873 0.00290 
AS24 0.022 0.029761 1.393 185295 0.00280 


0.0192 0.03248 | 0.000941 


0.0194 0.03411 | 0.001106 


AS25 0.024 0.027876 1.520 202140 0.00296 
AS26 0.018 0.032002 1.140 151605 0.00264 
AS31 0.020 0.03072 1.267 168450 0.00273 
0.008 AS32 0.021 0.032143 1.330 176873 0.00263 


0.0191 0.03151 | 0.000849 
0.0197 0.03723 | 0.001463 


0.0195 0.03544 | 0.001252 


0.0197 0.03743 | 0.001487 


AS33 0.023 0.032987 1.457 193718 0.00258 0.0198 0.03861 | 0.001637 


AS34 0.024 0.032344 1.520 202140 0.00262 0.0197 0.03771 | 0.001522 


Nikuradse Kb = 0.000080. 
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Table 9. Calculus of acrylic-sandpaper (Haland criteria), Scenario l. 


O de AS del El de 
AS1 | 0.017 | 0.027162 143 183 0.00315 0.0126 0.01812 0.03093 | 0.000621 
AS12 | 0.019 | 0.027302 160 028 0.00310 0.0127 0.01785 0.03126 | 0.000649 
0.001 | AS13 | 0.020 | 0.027369 168 450 0.00307 0.0127 0.01774 0.03141 | 0.000662 
AS14 | 0.022 | 0.027191 185 295 0.00306 0.0127 0.01753 0.03124 | 0.000653 
AS15 | 0.023 | 0.028838 1.457 193 718 | 0.00290 0.0129 0.01743 0.03367 | 0.000843 
AS21 | 0.018 | 0.0294330 1.140 151 605 0.00293 0.0129 0.01798 0.03428 | 0.000886 
AS22 | 0.020 | 0.0290041 1.267 168 450 0.00293 0.0129 0.01774 0.03377 | 0.000847 
0.004 | AS23 | 0.021 | 0.028582 1.330 176 873 0.00295 0.0128 0.01763 0.03321 | 0.000802 
AS24 | 0.022 | 0.029761 1.393 185 295 0.00284 0.0130 0.01753 0.03497 | 0.000952 
AS25 | 0.024 | 0.027876 1.520 202 140 0.00297 0.0128 0.01735 0.03232 | 0.000736 
AS26 | 0.018 | 0.032002 1.140 151 605 0.00274 0.0131 0.01798 0.03803 | 0.001232 
AS31 | 0.020 0.03072 1.267 168 450 0.00280 0.0130 0.01774 0.03627 | 0.001066 
0.008 | AS32 | 0.021 | 0.032143 1.330 176 873 0.00268 0.0131 0.01763 0.03841 | 0.001275 
AS33 | 0.023 | 0.032987 1.457 193 718 0.00260 0.0132 0.01743 0.03974 | 0.001416 
AS34 | 0.024 | 0.032344 1.520 202 140 0.00263 0.0132 0.01735 0.03884 | 0.001323 
Haland Kb = 0.000037. 
Table 10. Calculus of acrylic-sandpaper (Churchill criteria), Scenario 1. 
IE [Ea 1 O a le A 
ES AS11 0.017 | 0.027162 1.077 143 183 0.00291 0.0129 0.0167 0.03160 | 0.000740 
5 AS12 | 0.019 | 0.027302 1.203 160 028 0.00290 0.0129 0.0167 0.03180 | 0.000758 
E 0.001 | AS13 0.020 | 0.027369 1.267 168 450 0.00289 0.0129 0.0167 0.03190 | 0.000766 
e ASI4 | 0.022 | 0.027191 6895 185 295 0.00291 0.0129 0.0167 0.03164 | 0.000747 
S AS15 0.023 | 0.028838 1.457 193 718 0.00277 0.0130 0.0167 0.03403 | 0.000947 
E AS21 0.018 | 0.0294330 1.140 151 605 0.00273 0.0131 0.0167 0.03490 | 0.001021 
Bl AS22 | 0.020 | 0.0290041 1.267 168 450 0.00276 0.0130 0.0167 0.03428 | 0.000966 
E 0.004 | AS23 | 0.021 | 0.028582 1.330 176 873 0.00279 0.0130 0.0167 0.03366 | 0.000913 
E AS24 | 0.022 | 0.029761 1.393 185 295 0.00270 0.0131 0.0167 0.03538 | 0.001069 
E AS25 0.024 | 0.027876 1.5520 202 140 0.00285 0.0129 0.0167 0.03263 | 0.000828 
S AS26 | 0.018 | 0.032002 1.140 151 605 0.00255 0.0132 0.0167 0.03867 | 0.001396 
E AS31 0.020 0.03072 1.267 168 450 0.00263 0.0132 0.0167 0.03678 | 0.001202 
z 0.008 | AS32 | 0.021 | 0.032143 1.330 176 873 0.00254 0.0133 0.0167 0.03888 | 0.001420 
E AS33 0.023 | 0.032987 1.457 193 718 0.00249 0.0133 0.0167 0.04012 | 0.001559 
E AS34 | 0.024 | 0.032344 1.520 202 140 0.00253 0.0133 0.0167 0.03917 | 0.001454 
3, Churchill Kb = 0.000066. 
E 
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Table 11. Calculus of acrylic-sandpaper (Swamee criteria), Scenario l. 


S Code O (m/s) A e Re a dE Ab Aw A 
ASI1 0.017 0.027162 1.077 143 183 0.00253 0.0133 0.01416 0.03277 | 0.000788 
AS12 0.019 0.027302 1.203 160 028 0.00252 0.0133 0.01416 0.03298 | 0.000812 
0.001 AS13 0.020 0.027369 1.267 168 450 0.00252 0.0133 0.01416 0.03307 | 0.000822 
AS14 0.022 0.027191 1.393 185 295 0.00253 0.0133 0.01416 0.03281 | 0.000805 
AS15 0.023 0.028838 1.457 193 718 0.00241 0.0134 0.01416 0.03523 | 0.001020 
AS21 0.018 0.0294330 1.140 151 605 0.00237 0.0134 0.01416 0.03611 | 0.001091 
AS22 0.020 0.0290041 1.267 168 450 0.00240 0.0134 0.01416 0.03548 | 0.001037 
0.004 AS23 0.021 0.028582 1.330 176 873 0.00243 0.0134 0.01416 0.03486 | 0.000981 
AS24 0.022 0.029761 1.393 185 295 0.00235 0.0134 0.01416 0.03660 | 0.001149 
AS25 0.024 0.027876 1.520 202 140 0.00248 0.0133 0.01416 0.03382 | 0.000895 
AS26 0.018 0.032002 1.140 151 605 0.00222 0.0136 0.01416 0.03992 | 0.001488 
AS31 0.0200 0.03072 1.267 168 450 0.00229 0.0135 0.01416 0.03801 | 0.001288 
0.008 AS32 0.0210 0.032143 1.330 176 873 0.00221 0.0136 0.01416 0.04012 | 0.001520 
AS33 0.0230 0.032987 1.457 193 718 0.00216 0.0136 0.01416 0.04138 | 0.001668 
AS34 0.0240 0.032344 1.520 202 140 0.00220 0.0136 0.01416 0.04042 | 0.001560 


Swamee Kb = 0.000031. 


Scenario II. If geometry, absolute roughness 
materials (K,, K_) and influence perimeters (P,, 
P,,) are known, it is desired to calculate total 
roughness factor A of the model tunnel. 

The absolute roughness values K, , calculated 
in Scenario 1 (tables 8, 9, 10 and 11) are used in 
order to calculate A and the value of a Es with 
Eq. (17) (and Eq. 14). This predicted values of 
absolute roughness calculated with Colebrook's 
formula, let us to compare the measured values 
(tables 5, 6, 7 and the predicted values with the 
theoretical formulation). 

Table 12 show the results for Nikuradse, 
table 13 for Haland, table 14 for Churchill and 
table 15 for Swamee criteria, using the results 
obtained from scenario l. 


Statistical Comparison 


The statistical comparison of any friction factor 
equation with the Colebrook's equation can be 
done with the following procedure. 


+ Calculate the friction factor Roeii DY 
the criteria selected (Nikradse, Haland, 
Churchill, Swamee). 


Calculate the friction factor value K 


mC,i 


measured with the Colebrook's equation. 


Mean relative error 


a A 


10 
meanREK=- 


Maximal positive error 


maxRE*=max 


Ric _ PP. 
Rimci 


Maximal negative error 


maxRE” = max R 


Keredi E noi 
1 


Correlation ratio 


yx mc edi ] 


, 
e= 1-2 
” 2 


2 [Ko Ea) 


Calculate the following parameters: 


(Q1) 


(22) 


(23) 


(24) 
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Table 12. Calculus of acrylic-sandpaper (Nikuradse criteria), scenario II. 


oo la lao laa a ae a 
AS11 0.017 | 0.000762 | 1.077 | 143183 | 0.003248 | 0.0125 | 0.0187 | 0.0307 | 0.02726 | 0.0004616 
AS12 0.019 | 0.000778 | 1.203 | 160028 | 0.003237 | 0.0126 | 0.0187 | 0.0309 | 0.02740 | 0.0004700 
0.001 | AS13 0.020 | 0.000787 | 1.267 | 168450 |0.003232 | 0.0126 | 0.0187 | 0.0310 | 0.02747 | 0.0004741 
AS14 0.022 | 0.000765 | 1.393 | 185295 | 0.003246 | 0.0125 | 0.0187 | 0.0307 | 0.02729 | 0.0004634 
AS15 0.023 | 0.000976 | 1.457 | 193718 | 0.003118 | 0.0127 | 0.0189 | 0.0330 | 0.02894 | 0.0005674 
AS21 0.018 | 0.001059 | 1.140 | 151605 |0.003075 | 0.0127 | 0.0190 | 0.0339 | 0.02954 | 0.0006078 
AS22 0.020 | 0.000998 | 1.267 | 168450 |0.003106 | 0.0127 | 0.0189 | 0.0333 |0.029111 | 0.0005785 
0.004 | AS23 0.021 | 0.000941 | 1.330 | 176873 |0.003137 | 0.0127 | 0.0189 | 0.0327 | 0.02869 | 0.0005505 
AS24 0.022 | 0.001106 | 1.393 | 185295 |0.003052 | 0.0127 | 0.0190 | 0.0343 | 0.02987 | 0.0006309 
AS25 0.024 | 0.000849 | 1.520 | 202140 |0.003191 | 0.0126 | 0.0188 | 0.0317 | 0.027977 | 0.0005053 
AS26 0.018 | 0.001463 | 1.140 | 151605 |0.002906 | 0.0129 | 0.0192 | 0.0375 | 0.03213 | 0.0008005 
AS31 | 0.0200 | 0.001252 | 1.267 | 168450 |0.002988 | 0.0128 | 0.0191 | 0.0357 |0.030835 | 0.0007005 
0.008 | AS32 | 0.0210 | 0.001487 | 1.330 | 176873 |0.002898 | 0.0129 | 0.0192 | 0.0377 | 0.03227 | 0.0008119 
AS33 | 0.0230 | 0.001637 | 1.457 | 193718 |0.002848 | 0.0129 | 0.0193 | 0.0389 | 0.03312 | 0.0008820 
AS34 | 0.0240 | 0.001522 | 1.520 | 202140 |0.002886 | 0.0129 | 0.0193 | 0.0380 |0.032475 | 0.0008284 
Nikuradse Kb = 0.000080. 
Table 13. Calculus of acrylic-sandpaper (Haland criteria), scenario II. 
A A A A E E 5 
SS ASi1 0.017 |0.000621 | 1.077 | 143183 | 0.003222 | 0.0126 | 0.01812 | 0.03042 | 0.02688 | 0.000391 
a AS12 | 0.019 | 0.000649 | 1.203 | 160028 | 0.003169 | 0.0126 | 0.01785 | 0.03073 | 0.02702 | 0.000403 
E 0.001 | AS13 0.020 | 0.000662 | 1.267 | 168450 | 0.003145 | 0.0126 | 0.01774 | 0.03088 | 0.02708 | 0.000409 
5 ASI14 | 0.022 | 0.000653 | 1.393 | 185295 | 0.003130 | 0.0127 | 0.01753 | 0.03072 | 0.02691 | 0.000403 
E AS15 0.023 | 0.000843 | 1.457 | 193718 |0.002972 | 0.0128 | 0.01743 | 0.03310 | 0.02853 | 0.000504 
E AS21 0.018 | 0.000886 | 1.140 | 151605 | 0.003002 | 0.0128 | 0.01798 | 0.03370 | 0.02912 | 0.000533 
E AS22 | 0.020 | 0.000847 | 1.267 | 168450 | 0.003002 | 0.0128 | 0.01774 | 0.03320 | 0.028696 | 0.000509 
E 0.004 | AS23 0.021 | 0.000802 | 1.330 | 176873 | 0.003022 | 0.0128 | 0.01763 | 0.03265 | 0.02828 | 0.000484 
E AS24 | 0.022 | 0.000952 | 1.393 | 185295 | 0.002911 | 0.0129 | 0.01753 | 0.03437 | 0.02944 | 0.000563 
E AS25 0.024 | 0.000736 | 1.520 | 202140 | 0.003041 | 0.0127 | 0.01735 | 0.03177 | 0.027581 | 0.000445 
5 AS26 | 0.018 | 0.001232 | 1.140 | 151605 | 0.002808 | 0.0130 | 0.01798 | 0.03738 | 0.03165 | 0.000717 
E AS31 | 0.0200 | 0.001066 | 1.267 | 168450 | 0.002867 | 0.0129 | 0.01774 | 0.03564 | 0.030385 | 0.000626 
E 0.008 | AS32 | 0.0210 | 0.001275 | 1.330 | 176873 | 0.002751 | 0.0130 | 0.01763 | 0.03774 | 0.03179 | 0.000734 
E AS33 | 0.0230 | 0.001416 | 1.457 | 193718 | 0.002669 | 0.0131 | 0.01743 | 0.03905 | 0.03262 | 0.000804 
E AS34 | 0.0240 | 0.001323 | 1.520 | 202140 | 0.002700 | 0.0131 | 0.01735 | 0.03816 | 0.031988 | 0.000754 
3, Haland Kb = 0.000037. 
E 
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Table 14. Calculus of acrylic-sandpaper (Churchill criteria), scenario II. 


oO IE CA O O IE ES 0 
ASI1 0.017 [|0.000740 | 1.077 | 143183 |0.002968 | 0.0128 | 0.01666 | 0.031355 | 0.02666 | 0.000408 
AS12 | 0.019 |0.000758 | 1.203 | 160028 |0.002954 | 0.0128 | 0.01666 | 0.031584 | 0.02686 | 0.000422 
0.001 | AS13 | 0.020 |0.000766 | 1.267 | 168450 | 0.002947 | 0.0128 | 0.01666 | 0.031693 | 0.02696 | 0.000428 
ASI4 | 0.022 |0.000747 | 1.393 | 185295 |0.002962 | 0.0128 | 0.01666 | 0.031449 | 0.02681 | 0.000421 
AS15 | 0.023 |0.000947 | 1.457 | 193718 | 0.002822 | 0.0130 | 0.01666 | 0.033872 | 0.02852 | 0.000526 
AS21 | 0.018 |0.001021| 1.140 | 151605 | 0.002777 | 0.0130 | 0.01666 | 0.034708 | 0.02909 | 0.000559 
AS22 | 0.020 |0.000966 | 1.267 | 168450 | 0.002810 | 0.0130 | 0.01666 | 0.034094 | 0.028670 | 0.000533 
0.004 | AS23 | 0.021 |0.000913| 1.330 | 176873 |0.002843 | 0.0129 | 0.01666 | 0.033482 | 0.02825 | 0.000507 
AS24 | 0.022 |0.001069 | 1.393 | 185295 | 0.002751 | 0.0130 | 0.01666 | 0.035224 | 0.02948 | 0.000589 
AS25 | 0.024 |0.000828| 1.520 | 202140 | 0.002901 | 0.0129 | 0.01666 | 0.032467 | 0.027563 | 0.000466 
AS26 | 0.018 |0.001396 | 1.140 | 151605 | 0.002594 | 0.0132 | 0.01666 | 0.038502 | 0.03171 | 0.000750 
AS31 | 0.0200 | 0.001202 | 1.267 | 168450 | 0.002681 | 0.0131 | 0.01666 | 0.036617 | 0.030426 | 0.000655 
0.008 | AS32 | 0.0210 | 0.001420| 1.330 | 176873 | 0.002584 | 0.0132 | 0.01666 | 0.038731 | 0.03191 | 0.000768 
AS33 | 0.0230 | 0.001559| 1.457 | 193718 | 0.002531 | 0.0133 | 0.01666 | 0.039992 | 0.03280 | 0.000841 
AS34 | 0.0240 | 0.001454| 1.520 | 202140 | 0.002571 | 0.0134 | 0.01666 | 0.039045 | 0.032152 | 0.000789 
Churchill Kb = 0.000066. 
Table 15. Calculus of acrylic-sandpaper (Swamee criteria), scenario IL. 
E A O 0 
ASI1 0.017 [0.000788 | 1.077 | 143183 | 0.002607 | 0.0132 | 0.01416 | 0.03198 | 0.02671 | 0.000366 
AS12 | 0.019 [|0.000812 | 1.203 | 160028 | 0.002591 | 0.0132 | 0.01416 | 0.03227 | 0.02691 | 0.000383 
0.001 | AS13 | 0.020 |0.000822 | 1.267 | 168450 | 0.002584 | 0.0132 | 0.01416 | 0.03240 | 0.02700 | 0.000391 
ASI14 | 0.022 |0.000805 | 1.393 | 185295 | 0.002595 | 0.0132 | 0.01416 | 0.03219 | 0.02689 | 0.000389 
AS15 | 0.023 [|0.001020 | 1.457 | 193718 | 0.002468 | 0.0133 | 0.01416 | 0.03470 | 0.02854 | 0.000492 
AS21 | 0.018 |0.001091 | 1.140 | 151605 | 0.002432 | 0.0134 | 0.01416 | 0.03547 | 0.02911 | 0.000518 
AS22 | 0.020 |0.001037 | 1.267 | 168450 | 0.002459 | 0.0133 | 0.01416 | 0.03488 | 0.028695 | 0.000496 
0.004 | AS23 | 0.021 |0.000981 | 1.330 | 176873 | 0.002489 | 0.0133 | 0.01416 | 0.03427 | 0.02828 | 0.000471 
AS24 | 0.022 |0.001149 | 1.393 | 185295 | 0.002404 | 0.0134 | 0.01416 | 0.03607 | 0.02951 | 0.000555 
AS25 | 0.024 [|0.000895| 1.520 | 202140 | 0.002538 | 0.0133 | 0.01416 | 0.03327 | 0.027583 | 0.000434 
AS26 | 0.018 |0.001488 | 1.140 | 151605 | 0.002267 | 0.0135 | 0.01416 | 0.03936 | 0.03174 | 0.000709 
AS31 | 0.0200 | 0.001288 | 1.267 | 168450 | 0.002344 | 0.0134 | 0.01416 | 0.03747 | 0.030458 | 0.000618 
0.008 | AS32 | 0.0210 | 0.001520| 1.330 | 176873 | 0.002257 | 0.0135 | 0.01416 | 0.03965 | 0.03193 | 0.000732 
AS33 | 0.0230 | 0.001668 | 1.457 | 193718 | 0.002208 | 0.0136 | 0.01416 | 0.04096 | 0.03282 | 0.000808 
AS34 | 0.0240 | 0.001560| 1.520 | 202140 | 0.002243 | 0.0135 | 0.01416 | 0.04001 | 0.032176 | 0.000757 


Swamee Kb = 0.000031. 
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Standard deviation 


mCji pred,i 
i=1 K 
mC)ji 
Asp = 25 
SD a (25) 
Mean variation 
Koi 

Med _ m pate (26) 

Kim 


Where K_ is the average value of K' .. for the 
complete set of values: 


K,- Es) (27) 


In the table 16 is made a comparison with the 
mean relative, maximal positive, maximal nega- 
tive and mean variation error for each criteria 
(Nikuradse, Haland, Churchill, and Swamee). 


The maximal positive relative error maxRE* 
= 0.0936 is for Swamee criteria; the maximum 
negative error belongs to the Nikuradse criteria 
maxRE' = 0.144. 

In table 20 is shown a summary of this cal- 
culus; in which are shown the correlation ratio, 
standard deviation and mean variation of the 
error. 

The statistical analysis, show that the best 
correlation ratio and the minimum standard 
deviation is for Churchill criteria. 

The minimum difference in the in the mean 


.. _ Rod 


value A,, = is for Haland criteria 


(-0.03665) although the obtained for Churchill 
(0.05517) and Swamee (0.05515) are very similar. 


Conclusions 


It was proposed the Eq. (8) thatis a Colebrook- 
White form, and let to determine the absolute 
roughness in this kind of tunnels geometries; 
even when it is applied to the composite 
roughness. In addition, it was reviewed other 


Table 16. Calculus of acrylic-sandpaper (Nikuradse criteria). 


2 
da men Kmeas Knik a maxRE* | maxRE (K, -K. y (K, - K, A Bo ) 

0.0010 | 0.0170 1.077 0.000404 | 0.000462 0.144 -0.144 0.144 | 3.363E-09 | 1.629E-07 | 2.064E-02 
0.0190 1.203 0.000418 | 0.000470 0.126 -0.126 0.126 | 2.754E-09 | 1.743E-07 | 1.580E-02 
0.0200 1.267 0.000424 | 0.000474 0.118 -0.118 0.118 2.515E-09 | 1.797E-07 | 1.399E-02 
0.0220 1.393 0.000417 | 0.000463 0.111 -0.111 0.111 2.127E-09 | 1.741E-07 | 1.222E-02 
0.0230 1.457 0.000523 | 0.000567 0.085 -0.085 0.085 1.988E-09 | 2.734E-07 | 7.273E-03 

0.0040 | 0.0180 1.140 0.000553 | 0.000608 0.100 -0.100 0.100 | 3.027E-09 | 3.056E-07 | 9.904E-03 
0.0200 1.267 0.000528 | 0.000579 0.095 -0.095 0.095 | 2.522E-09 | 2.791E-07 | 9.036E-03 
0.0210 1.330 0.000502 | 0.000550 0.096 -0.096 0.096 | 2.315E-09 | 2.524E-07 | 9.174E-03 
0.0220 1.393 0.000584 | 0.000631 0.080 -0.080 0.080 |2.163E-09 | 3.416E-07 | 6.333E-03 
0.0240 1.520 0.000462 | 0.000505 0.093 -0.093 0.093 1.835E-09 | 2.138E-07 | 8.583E-03 

0.0080 | 0.0180 1.140 0.000745 | 0.000800 0.075 -0.075 0.075 3.094E-09 | 5.548E-07 | 5.576E-03 
0.0200 1.267 0.000650 | 0.000701 0.078 -0.078 0.078  |2.561E-09 | 4.224E-07 | 6.064E-03 
0.0210 1.330 0.000763 | 0.000812 0.065 -0.065 0.065 2.424FE-09 | 5.817E-07 | 4.167E-03 
0.0230 1.457 0.000836 | 0.000882 0.055 -0.055 0.055  |2.150E-09 | 6.983E-07 | 3.078E-03 
0.0240 1.520 0.000784 | 0.000828 0.057 -0.057 0.057 1.981E-09 | 6.145E-07 | 3.225E-03 

Mean 0.000573 | 0.000622 0.092 -0.055 0.144 

bd 0.99647 
Asp 0.009004 
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Table 17. Calculus of acrylic-sandpaper (Haland criteria). 


p 
Ss o es Kmeas Khal pe e maxRE* | maxRE ( K,,- nor (K,, - K,) (E) 
0.0010 | 0.0170 1.077 | 0.000404 | 0.000391 0.03230 0.03230 | -0.03230 | 1.700E-10 | 1.629E-07 | 2.074E-02 
0.0190 1.203 | 0.000418 | 0.000403 0.03388 0.03388 | -0.03388 | 2.001E-10 | 1.743E-07 | 1.148E-03 
0.0200 1.267 | 0.000424 | 0.000409 0.03452 0.03452 | -0.03452 | 2.142E-10 | 1.797E-07 | 1.192E-03 
0.0220 1.393 | 0.000417 | 0.000403 0.03527 0.03527 | -0.03527 | 2.167E-10 | 1.741E-07 | 1.244E-03 
0.0230 1.457 | 0.000523 | 0.000504 0.03700 0.03700 | -0.03700 | 3.743E-10 | 2.734E-07 | 1.369E-03 
0.0040 | 0.0180 1.140 | 0.000553 | 0.000533 0.03561 0.03561 | -0.03561 | 3.876E-10 | 3.056E-07 | 1.268E-03 
0.0200 1.267 | 0.000528 | 0.000509 0.03611 0.03611 -0.03611 | 3.639E-10 | 2.791E-07 | 1.304E-03 
0.0210 1.330 | 0.000502 | 0.000484 0.03617 0.03617 | -0.03617 | 3.302E-10 | 2.524E-07 | 1.308E-03 
0.0220 1.393 | 0.000584 | 0.000563 0.03717 0.03717 | -0.03717 | 4.719E-10 | 3.416E-07 | 1.382E-03 
0.0240 1.520 | 0.000462 | 0.000445 0.03668 0.03668 | -0.03668 | 2.877E-10 | 2.138E-07 | 1.345E-03 
0.0080 | 0.0180 1.140 | 0.000745 | 0.000717 0.03684 0.03684 | -0.03684 | 7.529E-10 | 5.548E-07 | 1.357E-03 
0.0200 1.267 | 0.000650 | 0.000626 0.03699 0.03699 | -0.03699 | 5.780E-10 | 4.224E-07 | 1.369E-03 
0.0210 1.330 | 0.000763 | 0.000734 0.03755 0.03755 | -0.03755 | 8.203E-10 | 5.817E-07 | 1.410E-03 
0.0230 1.457 | 0.000836 | 0.000804 0.03793 0.03793 | -0.03793 | 1.005E-09 | 6.983E-07 | 1.439E-03 
0.0240 1.520 | 0.000784 | 0.000754 0.03807 0.03807 | -0.03807 | 8.906E-10 | 6.145E-07 | 1.449E-03 
Mean 0.000573 | 0.000552 0.03614 0.03807 | -0.03230 
e 0.999324 
Asp 0.00262 
Table 18. Calculus of acrylic-sandpaper (Churchill criteria). 
2 
S a MA Kmeas Kchur Ea maxRE* | maxRE' | (K,-K, ñ (K, - K,)' E] 
(m/s) | (m/s) Reni KE 
0.0010 | 0.0170 | 1.077 | 0.000404 0.000408 0.01181 -0.01181 0.01181 | 2.272E-11 | 1.629E-07 | 1.395E-04 
0.0190 | 1.203 | 0.000418 0.000422 0.01032 -0.01032 | 0.01032 | 1.857E-11 | 1.743E-07 | 1.065E-04 
0.0200 | 1.267 | 0.000424 0.000428 0.00976 -0.00976 | 0.00976 | 1.712E-11 | 1.797E-07 | 9.526E-05 
0.0220 | 1.393 | 0.000417 0.000421 0.00918 -0.00918 | 0.00918 | 1.468E-11 | 1.741E-07 | 8.430E-05 
0.0230 | 1.457 | 0.000523 0.000526 0.00526 -0.00526 | 0.00526 | 7.565E-12 | 2.734E-07 | 2.767E-05 
0.0040 | 0.0180 | 1.140 | 0.000553 0.000559 0.01185 -0.01185 0.01185 | 4.289E-11 | 3.056E-07 | 1.403E-04 
0.0200 | 1.267 | 0.000528 0.000533 0.00978 -0.00978 | 0.00978 | 2.667E-11 | 2.791E-07 | 9.556E-05 
0.0210 | 1.330 | 0.000502 0.000507 0.01009 -0.01009 | 0.01009 | 2.570E-11 | 2.524E-07 | 1.019E-04 
0.0220 | 1.393 | 0.000584 0.000589 0.00735 -0.00735 | 0.00735 | 1.845E-11 | 3.416E-07 | 5.402E-05 
0.0240 | 1.520 | 0.000462 | 0.000466 0.00857 -0.00857 | 0.00857 | 1.569E-11 | 2.138E-07 | 7.339E-05 
0.0080 | 0.0180 | 1.140 | 0.000745 0.000750 0.00674 -0.00674 | 0.00674 | 2.518E-11 | 5.548E-07 | 4.538E-05 
0.0200 | 1.267 | 0.000650 0.000655 0.00712 -0.00712 | 0.00712 | 2.144E-11 | 4.224E-07 | 5.076E-05 
0.0210 | 1.330 | 0.000763 0.000768 0.00646 -0.00646 | 0.00646 | 2.428E-11 | 5.817E-07 | 4.175E-05 
0.0230 | 1.457 | 0.000836 | 0.000841 0.00622 -0.00622 | 0.00622 | 2.699E-11 | 6.983E-07 | 3.864E-05 
0.0240 | 1.520 | 0.000784 0.000789 0.00635 -0.00635 | 0.00635 | 2.479E-11 | 6.145E-07 | 4.034E-05 
Mean 0.000573 | 0.00054139 0.00846 -0.00526 | 0.01185 
e 0.9999682 
A 7.5682E-05 
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Table 19. Calculus of acrylic-sandpaper (Swamee criteria). 


2 
Q V Kouci ES Road ES ER 2 2 K,, = K. 
(m/s) A Kmeas Kswam RE maxRE* | maxRE (K, - K.) (K, a K.,) K 


mCji 


0.0010 | 0.0170 | 1.077 | 0.000404 | 0.000366 0.0936 0.0936 -0.0936 
0.0190 | 1.203 | 0.000418 | 0.000383 0.0818 0.0818 -0.0818 
0.0200 | 1.267 | 0.000424 | 0.000391 0.0769 0.0769 -0.0769 
0.0220 | 1.393 | 0.000417 | 0.000389 0.0669 0.0669 -0.0669 
0.0230 | 1.457 | 0.000523 | 0.000492 0.0582 0.0582 -0.0582 
0.0040 | 0.0180 | 1.140 | 0.000553 | 0.000518 0.0631 0.0631 -0.0631 
0.0200 | 1.267 | 0.000528 | 0.000496 0.0616 0.0616 -0.0616 
0.0210 | 1.330 | 0.000502 | 0.000471 0.0619 0.0619 -0.0619 
0.0220 | 1.393 | 0.000584 | 0.000555 0.0506 0.0506 -0.0506 
0.0240 | 1.520 | 0.000462 | 0.000434 0.0609 0.0609 -0.0609 
0.0080 | 0.0180 | 1.140 | 0.000745 | 0.000709 0.0480 0.0480 -0.0480 
0.0200 | 1.267 | 0.000650 | 0.000618 0.0498 0.0498 -0.0498 
0.0210 | 1.330 | 0.000763 | 0.000732 0.0400 0.0400 -0.0400 
0.0230 | 1.457 | 0.000836 | 0.000808 0.0329 0.0329 -0.0329 
0.0240 | 1.520 | 0.000784 | 0.000757 0.0337 0.0337 -0.0337 


1.426E-09 | 2.864E-08 8.75E-03 
1.168E-09 | 2.412E-08 6.70E-03 
1.063E-09 | 2.217E-08 | 5.92E-03 
7.796E-10 | 2.419E-08 | 4.48E-03 
9.255E-10 | 2.498E-09 3.39E-03 
1.217E-09 | 4.004E-10 3.98E-03 
1.059E-09 | 1.983E-09 3.80E-03 
9.667E-10 | 4.968E-09 3.83E-03 
8.759E-10 | 1.343E-10 2.56E-03 
7.932E-10 | 1.219E-08 3.71E-03 
1.277E-09 | 2.960E-08 2.30E-03 
1.048E-09 | 5.940E-09 2.48E-03 
9.314E-10 | 3.604E-08 1.60E-03 
7.540E-10 | 6.907E-08 1.08E-03 
6.984E-10 | 4.454E-08 1.14E-03 


Mean 0.000573 | 0.0005414 0.0587 0.0936 -0.0329 
9 0.975248717 
Asp 0.003714213 
Table 20. Summary of statistical analysis. 
Nikuradse Haland Churchill Swamee 
o) Asp Ani e Asp Ani 0 Asp Aj o) Asp An 
0.99647 0.0090 -0.08551 | 0.99932 | 0.00262 | -0.03665 | 0.99997 |7.5682E-05 | 0.05515 |0.9752487 | 0.003714 | 0.05517 


equations like Nikuradse, Haland, Churchill 
and Swamee. 


In future research, will be possible to ana- 
lyze the hydraulic behavior in tunnels having 
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Elfman's (1993) and Marengo-Mogollón 
(2005) theoretical analysis was reviewed and 
was proposed two Scenario's analysis. They 
were compare with models measurements and 
they were validate. When comparing the results 
obtained from simple equations with those de- 
termined from the experiments their validity 
is demonstrated to calculate coefficients A, A, 
and K, (Scenario 1), and A, A, A and K,,. in 
Scenario II. The theoretical development is quite 
acceptable. 

Theoretical model can be applied in model 
development with this kind of analysis. 
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other geometries (horseshoe and circular cross 
sections with composite roughness), so as to be 
able to define the behavior of the tunnels when 
functioning as open channels. 

It will be possible to validate the equations 
proposed for the composite roughness with 
their respective geometries, and to investigate 
the behavior of the equal-velocity curves with 
tunnels operating as full pipes. 

It should be made special mention in order 
to measure in prototype tunnels and com- 
pare with the results obtained in hydraulic 

models. 
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With respect to the statistical analysis, stan- 
dard deviation (A,,), correlation ratio (9), and 
maximal relative errors (maxRE”), are quite low. 

The minimum values of standard deviation 
A; = 7.58 belongs to Churchill criteria besides 
the best correlation ratio O = 0.9999682. The 
maximal relative error belongs to Swamee 
maxRE* = 0.0936 and the maximum negative 
error belongs to the Nikuradse criteria maxRE” 
=0.144. 


The minimum difference in the mean value 


K_o-K 
, = 2 — PL is for Haland criteria (-0.03665) 


ao E obtained for Churchill (0.05517) 
and Swamee (0.05515) are very similar. 

Under this analysis it is possible to say that 
the best theoretical criteria in order to obtain 
an accurate behavior is with Churchill criteria 
that offers the best correlation ratio, followed by 
Haland, Nikuradse and Swamee. 

If it is selected standard, division also is 
Churchill criteria. If mean criteria is selected, 
Haland gives the best approximation. 

The results obtained with Haland criteria are 
very good also. 


Notation 


A 
a 


cross sectional area. 


height of triangle defining the area of 

influence of the floor material in tunnels 

with composite roughness. 

A, and A,,= areas of influence of the material 
for the floor slab and for walls and vault, 
respectively. 

d = depth of water at tunnel entrance. 

D = pipe diameter. 

d/D = ratio between inlet hydraulic head with 
respect to the tunnel equivalent diameter. 

D, = equivalent diameter to the circular cross 
section. 

f = Marchi's shape factor. 

F  = total shear force. 

F, and F,,= Shear forces at floor, and at walls and 

vault, respectively. 

acceleration of gravity. 


=— 0 
1! 


lost of energy. 


K = absolute roughness of the material of the 
circular conduit. 

K, and K,, = roughness coefficient at floor and at 
walls and vault, respectively, in a tunnel 
with composite roughness. 


K, = conversion factor of units for Manning/'s 
formula. 
K, = equivalent absolute roughness of the 


material for non-circular conduits. 
1 = length of the section. 
L, = scale of lines. 
LP = area of section under study. 
n, = scale of roughness. 
P = perimeter. 
Q = flow rate. 
OQ = scale of flow rates. 
Re = Reynolds' number. 
Rh = hydraulic radius. 
R, = Reynolds” number calculated with the 
hydraulic radius. 
S = hydraulic gradient. 
S, and S,= slopes of tunnels. 
V = mean velocity. 
V, and V,, = mean values at floor, walls and 
vault. 
= scale of velocities. 
= maximum flow velocity. 
=real velocity registered with a Prndtl- 
Pitot tube. 
= density of fluid. 


p 

A = increment as a percentage. 

v = kinematic viscosity of water. 

T = Shear stress at the wall. 

A, n, C = coefficients of resistance to flow by 
Darcy-Weisbach, Manning and Chezy, 
respectively. 

1, and 1,, = shear stresses at perimeter of floor 
slab and at perimeter of walls and vault, 
respectively. 

A, and A, = coefficients of resistance to flow at 
floor, and at walls and vault, respectively. 

A. = coefficient of resistance calculated with 
programs developed by Marengo- 
Mogollón. 

Ah, and Ah,= hydraulic head losses at floor and 
at walls and vault, respectively; 

A, = average coefficient of resistance. 


m 
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Annex I. Roughness criteria for vault section. 


Table I. Values of the dimensionless coefficient A and absolute roughness to estimate head losses in 
pressurized conduits. 


Author Value ofA Value of K Eq. 
1 loe 3.71D Pipes K= 2A/P 
NN K 10 (0 -1,74) 2A/P 
Nikuradse A 111.11.11 
(1983) E ea Al 
qa 174+ og A Other shape 
%A KP 
A q il 0.9009 E 
Haaland (1983) 2 
(Franzini ée E k= (ES te > AZ 
Finnemore, 1999) SA Eon En (22) 
penal ci v 
V 
)= pl 
2 
Churchill-Barr K= (E) 107420 _ EA 
(Bombardelli € alo 2K , 5,76 2 7) A3 
García, 2003) 814,83 q [2voy” pi 
V 
a 0.25 
: A : 
Swamee-Jain ren 14,83b 19120 5.76 
(1976) 2K 5,76 a 2 5 py e 
(Bombardelli é2 ans (22) 
García, 2003) ó e 
Y 
Colebrook 
1 2K. 252 14.830). | 2h 2.52 
(Franzini de yn = -2log + Ko 22 119424 _ AS 
ya 14.83b 2Vb 2 2Vb 
Finnemore, 1999) Y A 
a 1/8 

¡SS Guo-Pierre 0.3164 R 
on a ITA sd 
(LAR, 2004) : as 
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Abacus A vs. Re Tunnels composite roughness 


A 
0.090 — e AS - e 
0.085 a 
0.080 “== TT TT] A 4er Acrylic-carpet HH 
PTTTTH z (experimenta) HO 0.06 
0.075 O 
Wo petytitrar pet (calculated) 
a 
0.070 > 
0.065 '.N LL 904 
0.060 Turbulent zqne, 
0.055 
0.050 dea 
0.01 
0.040 ¡MA ll 
rylic-plkstic (experimental) LEADA 
Lic de 0.006 
RL Acrylic-plpstic (calculated) || 
f=64/Re Y jo 
0.030 Poiseuille PE 0.004 
laminar flow ; A 
0.025 Se 
1000 10 000 100 000 1 000 000 10 000 000 100 000 000 


Re 


Figure 5. Abacus of As. R, curves for all tunnels with uniform and composite roughness. 
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